Integrated multi vector vortex beam generator 



Sebastian A. Schulz, Taras Machula, and Ebrahim Karim£| 

Department of Physics, University of Ottawa, 150 Louis Pasteur, Ottawa, Ontario, KIN 6N5 Canada 



CO 

o 

CN 



o 

H— > . 

Or 

q 
o 

Or 



> 

OO 

o 



x 



Robert W. Boyd 

Department of Physics, University of Ottawa, 150 Louis Pasteur, Ottawa, Ontario, KIN 6N5 Canada and 
Institute of Optics, University of Rochester, Rochester, New York, 14627, USA 

A novel method to generate and manipulate vector vortex beams in an integrated, ring resonator 
based geometry is proposed. We show numerically that a ring resonator, with an appropriate grating, 
addressed by a vertically displaced access waveguide emits a complex optical field. The emitted beam 
possesses a specific polarization topology, and consequently a transverse intensity profile and orbital 
angular momentum. We propose a combination of several concentric ring resonators, addressed 
with different bus guides, to generate arbitrary orbital angular momentum audit states, which could 
potentially be used for classical and quantum communications. Finally, we demonstrate numerically 
that this device works as an orbital angular momentum sorter with an average cross-talk of — lOdB 
between different orbital angular momentum channels. 



INTRODUCTION 

It is well known that a light beam carries a linear 
momentum, which can be transferred to an absorbing or 
reflecting medium [l| . This linear momentum of light is 
proportional to the time-average of the Poynting vector. 
The Poynting vector associated with a plane wave is 
purely longitudinal and has no transverse components, 
consequently the linear momentum is parallel to the 
propagation direction. However, since a plane wave 
carries infinite energy it is not physically realizable. 
Therefore, all generated optical beams are solutions of 
the paraxial wave equation, leading to a non-vanishing 
component of the linear momentum in the transverse 
plane 0. Such a paraxial beam might possess a net 
orbital angular momentum (OAM) in the propagation 
direction as a consequence of having a finite transverse 
linear momentum [l], Stationary azimuthal linear 
momentum leads to a quantized OAM of Ifi per photon 
in the propagation direction, where I and h are an integer 
and the reduced Planck constant respectively [1, The 
phase front associated with such a beam forms a helical 
shape of exp (iitp) , where ip is the azimuth angle of the 
polar coordinate system. Additionally, light possesses 
an independent form of angular momentum, called spin 
angular momentum (SAM). The SAM is related to 
the vectorial property of optical fields, with a value of 
±S per photon, corresponding to left and right hand 
circularly polarized light, respectively @, [tJ- Although 
the SAM and OAM are independent features of the 
optical field, they can interact with each other under 
certain circumstances 

Many methods to generate and manipulate optical OAM 
have been presented in the past, due to the potentially 
use of OAM in a wide variety of applications, such as 
optical lithography, astronomy, optofluidics and optical 
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telecommunication (llHl^. Additionally to the afore- 
mentioned classical applications, optical OAM opens 
up a promising perspective in quantum communication, 
where a audit state can be used to encode informa- 
tion [17H20| . Implementing a single photon qudit state 
instead of a multi photon quNit (N=d) state results in 
a significantly lower power consumption during state 
preparation, transmission and detection processes. The 
previously suggested approaches to generate and ma- 
nipulate the optical OAM include computer-generated 
holograms screened on a spatial light modulator (SLM), 
astigmatic mode converters, spiral phase plates and 
spin-to-orbit conversion in inhomogeneous birefringent 
plates (g-plates) 0, [UHlll ■ These techniques, however, 
have several limitations. For instance, beam generation 
using SLMs has a low efficiency and the beam quality 
is bounded by the pixel size of nematic liquid crystal 
cells in the spatial light modulator, while the other 
approaches are stationary and cannot be dynamically 
controlled. 

Integrated photonics for quantum computation has 
generated significant interest, as it allows for a re- 
duced footprint and energy consumption compared to 
bulk optics components, as well as improved stabil- 
ity [24ll25| . Additionally, silicon photonics based devices 
are complementary metal-oxide-semiconductor (CMOS) 
compatible, offering the opportunity of large scale 
fabrication and integration with electronic circuits. Y. 
F. Yu et al. proposed an integrated photonics model to 
generate optical OAM, usin g a ring resonator surrounded 
by a group of nano-rods [26|. The ring resonator is 
side coupled to a bus waveguide, with nano-rods at 
the outer circumference of the ring resonator acting as 
scatterers for the local evanescent field. Each nano-rod 
has a specific phase delay with respect to neighboring 
nano-rods, leading to a helical phase front and therefore 
a finite OAM. More recently, X. Cai et al. fabricated 
a compact micro optical vortex beam emitter on a 
silicon-on-insulator chip, where the nano-rods proposed 
by Y. F. Yu. et al. were replaced by an in-ring angular 
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Figure 1. Schematic of proposed configuration of a silicon 
bus wave-guide and ring resonator. The ring resonator is on 
the top of the access waveguide, separated by a 275 nm thick 
layer of silica. The angular grating inside the ring resonator 
are similar to that suggested in Ref. [27l |. 



grating [27] . In fact, the beam generated by either 
angular nano-rods or a grating does not carry a pure 
OAM value, instead the emitted light beam owns a 
space-variant polarization pattern. While providing a 
significant reduction in footprint when compared to 
other methods for the generation of OAM beams, these 
methods are not suitable for the generation of OAM 
qudit states, as a single ring resonator will only emit 
"vector vortex" beams with a bounded OAM value for a 
given wavelength. 

In this letter, we propose a novel design for the gener- 
ation of vector vortex beams, where the bus waveguide 
and the ring resonator are located in different planes (2|| . 
Numerical simulations show that this out of plane con- 
figuration leads to very good coupling between the bus 
guide and the ring resonator, comparable to the side- 
coupled configuration, while the quality of the generated 
vector vortex beam is not affected by the bus waveg- 
uide placement. Furthermore, the proposed design gives 
the opportunity to build up concentric ring resonators 
with different radii that can be addressed independently 
through individual access waveguides. We show numeri- 
cally that adjusting both the phase and the amplitude of 
the bus waveguides input can be used to encode a qudit 
state. 



THEORETICAL MODEL 

As discussed in the introduction, it has been shown 
that a sequence of nano-rods or notches on a waveguide 
can be used as a grating to scatter the evanescent wave 
near a ring resonator (2lj ■ Under the correct condi- 
tions the emitted beam will carry OAM, with the OAM 
value being determined by the difference between the 
resonator mode number and the number of scattering 
elements. In the previous configurations both the ac- 
cess waveguide and the ring resonator were in the same 
plane, with a horizontal (side) coupling between them. 
This horizontal coupling configuration limits the devices 
to mainly preliminary applications, since the higher di- 
mensions of the OAM Hilbert space stay inaccessible. In 
order to overcome this issue, we propose a configuration 
where the waveguide and ring resonator are in a sepa- 
rate plane, as shown in Fig. ((T|). In the proposed scheme 
the access waveguide is still coupled to the ring resonator 
through evanescent coupling, albeit in the vertical direc- 
tion. On resonance the light inside the ring resonator 
satisfies the equation 



R Q = m\/(2im e ff), 



(1) 



where m, n e ff and A stand for the ring resonator mode 
number (a positive integer), the effective index of the 
waveguide mode and the free-space wavelength, respec- 
tively. The ring resonator mode number, to, is denned 
as the number of optical periods in the resonator cir- 
cumference, 2ttRq. This standing wave possesses a "net" 
azimuthal linear momentum of /3 m = 27m e ///A = m/Ro 
in the plane of ring resonator. In order to couple this 
beam out of the ring resonator (into the vertical direc- 
tion) , we choose an appropriate grating with elements on 
the inside wall of the ring resonator and equal angular 
spacing, similar to those suggested in Ref. [27| . 
For the remainder of this work we assume that the res- 
onator is on a resonance, with mode number m. The 
grating elements can be designed to be at any relative 
phase position of the standing wave. For instance, let 
us assume a ring resonator with m — 1, i.e. only one 
optical period is contained within the ring. Four sym- 
metric objects inside the ring resonator with a relative 
angle of 7r/2 could be used to pick up a beam with a 
relative phase of ir/2. However, adding additional scat- 
tering elements improves the quality of the phase front 
as well as the OAM state purity, but does not affect the 
beam structure. Instead, a higher order mode of the ring 
resonator should be excited in order to change the beam 
phase-front. To clarify this point, let us consider a ring 
resonator with m = 4 and the four symmetric scattering 
objects from the previous example. In this case, the scat- 
tering elements pick up the coupled wave at four different 
positions, all of which are in phase. Therefore, in the ab- 
sence of diffraction effects, the scattered beam possesses 
a uniform phase-font and consequently no OAM. More 
generally, it has been shown theoretically that the emit- 
ted wave from the grating interferes constructively if the 
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Figure 2. Theoretical intensity and polarization patterns for a beam emitted from a ring resonator. The upper row is for light 
polarized perpendicular to the waveguide, yet in the plane of the resonator, while the lower row is for light polarized parallel 
to the waveguide, leading to a radial and an azimuthal polarization distribution inside the ring resonator respectively. Each 
column corresponds to a specific angular phase matching condition and consequently value of OAM. Interestingly, for the case 
where £ = ±1, the intensity pattern of the outgoing beam does not form a doughnut shape, and its polarization switches from 
radial to azimuthal and vice versa. This type of beam is know as the Poincare or polarization-singular beam. 



following angular phase matching condition is satisfied; 

I = m - q, (2) 

where q stands for the number of grating elements (e.g., 
see supplementary material of Ref. \2T\). 
In addition to the phase property of light, the light in- 
side the ring resonator also has a well defined polariza- 
tion. The electric field can either lie in the same hori- 
zontal plane as the ring resonator, i.e. transverse electric 
(TE), or orthogonal to this plane, i.e. transverse mag- 
netic (TM) polarization. Note that these definitions of 
polarization are standard in integrated optics, however 
they differ from the standard optics definitions. In this 
work we exclusively use the integrated optics definitions, 
as given above. 

We assume that the grating is made of a linear, isotropic, 
and nonmagnetic medium. Under this circumstance and 
assuming a weak scattering regime the far-field projection 
(at r) of the optical fields (E, B) of the light scattered by 
each grating pitch (at ro = Ror) is given by: 

ikr 

E(r) - A(r ) 

r 

B(r) = - x A(r ) — , (3) 

c r 

where c and k = kk are the speed of light in vac- 
uum and the scattered beam wave-vector respectively. 

A(r ) oc — xk x kx E(r ) ) is the scattering ampli- 



tude, where \ is the grating susceptibility and E(ro) is 
the evanescent wave at the grating (29|. As we expected, 
both E and B obey the Maxwell equations and there- 
fore E, B and k are mutually orthogonal vectors. From 
this follows intuitively that for the case of a TM mode 
the light scattered by the grating is in the ring resonator 
plane and will not contribute to the out of plane radia- 
tion (in fact, the scattering amplitude A(ro) vanishes). 
Therefore, for the remainder of this work we will consider 
the case of TE polarization only, where the electric field 
inside the ring resonator as well as the evanescent wave 
are both either radially or azimuthally polarized, cor- 
responding to electric fields orthogonal to or along the 
propagation direction in the access waveguide. Hence, 
the scattered optical field of Eq. ^ is either radially or 
azimuthally polarized as well, i.e. A(ro) oc A(ro)r or 
A(ro) oc A(ro)</?, where r and (p are the radial and az- 
imuthal unit vectors in the polar coordinates. 
Therefore the overall scattered optical field has a complex 
structure in the far-field. Assuming the angular condi- 
tion stated in Eq.@, its phase-front forms a helix, de- 
scribed by e v , and it is either radially or azimuthally 
polarized, depending on the polarization of light in the 
access waveguide. Thus, we can express the optical field 
emerging form the ring resonator in one of the following 
forms: 

• radially polarized beam 
E r (r) =E (r,z)e u *r 
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Figure 3. Sketch showing the stereoscopic mapping of the 
polarization Poincare sphere on a beam transverse plane for 
the case of radially polarized beam with I = — 1. The beam 
possesses a right-handed circular polarization at the center 
(named C-point), tangent to the sphere's south pole. From 
this point the polarization transitions through right handed 
elliptical to linear polarization at the so called L-line, a map- 
ping of the Poincare sphere's equator. The upper hemisphere 
of the Poincare sphere, however, is mapped to the region out- 
side of the L-line, which therefore has a left handed polariza- 
tion (initially elliptical, then circular). In the case of a radi- 
ally polarized beam with I = +1, the stereoscopic mapping is 
on the north pole, consequently the polarization handedness 
changes from right to left handed. 



E (r,z) 
V2 



[\L) n \£-l) o +\R)Jt + l) ], (4) 



• azimuthally polarized beam 

E <p (r) = E (r,z)e i ^cp 
iE (r,z) 



V2 



\L)„\l-l) + \R)Ji+l) o },(5) 



More interestingly, under specific angular phase match- 
ing conditions, where £ = ±1, the emerging beam from 
the ring resonator forms a complex polarization structure 
with a non-zero intensity at the beam center, for both 
types of input polarization. These beams are known as 
a family of Poincare or polarization-singular beams f3l| — 
|33| . In fact, there is a one-to-one stereoscopic mapping 
between the polarization Poincare sphere and the beam 
polarization pattern on the transverse plane. However, 
for the case of £ = — 1, there is a stereoscopic mapping 
between a pole-to-pole path on the polarization Poincare 
sphere and the beam's transverse plane, as shown in Fig. 
([3]) . At the centre of the beam is a right-handed circular 
polarization. At a certain region, where the two ampli- 
tude of Eq. ((J]) are equal, the polarization has changed 
to a linear polarization and after this radius, the po- 
larization handedness changes into left-handed circular 
polarization. The Stokes parameters corresponding to 
Poincare beams change across the transverse beam plane. 
In certain regions or points, where the polarization is ei- 
ther a circular (Si and 52 = 0) or a linear (S3 = 0) 
polarization, the orientation of polarization ellipses or 
polarization handedness will be undefined. These points 
or region are called C-points or L-lines, respectively. In 
this analysis, we have neglected any possible variations in 
the radial beam profile due to confinement of the scatter- 
ing sources and instead a Laguerre-Gauss mode of radial 
order zero was assumed. Nonetheless, any paraxial beam 
as well as the emitted beam from the ring resonator can 
be expanded in terms of the Laguerre-Gauss modes, since 
they represent a complete set of solution to the paraxial 
wave equation. Therefore, and without loss of generality, 
we choose to truncate this expansion after the first or- 
der and have neglected the higher excited radial modes. 
Moreover, the polarization distribution pattern is altered 
dramatically by free-space prop agation, as has been ex- 
perimentally verified in Ref [33j. 



SIMULATION 



where Eq(t, z) is the radial amplitude of the scalar optical 
field, \S) V and \0) Q are the states of SAM and OAM in 
the Dirac notation, and L and R stand for left-circularly 
and right-circularly polarized beams respectively. The 
polarization and intensity patterns of the output beams 
for different angular phase matching conditions, i.e. i, 
and different polarization coupling conditions are shown 
in Fig. @. When the angular phase matching is set at 
£ = 0, the emitted beam from the ring resonator does not 
possess a net OAM value. However, it has either radial 
or azimuthal polarization, depending on the polarization 
of light in the access waveguide. These polarization con- 
ditions have a singularity at the centre and therefore the 
beam shape is a doughnut, as can be seen in the middle 
column of Fig. @. These optical beams have interest- 
ing applications in lithography and data storage, where 
they provide minimum spot size under tight focusing [30j | . 



In the previous section, we have discussed the proper- 
ties of a beam emitted from a single ring resonator with 
scattering elements. The discussion was completely gen- 
eral and no assumption about the geometrical arrange- 
ment of the ring resonator and access waveguide was 
made. However, as stated in the introduction, a single 
ring resonator cannot be used to generate a qudit state, 
and we instead propose a system of multiple concentric 
ring resonators that are addressed through vertically dis- 
placed access waveguides. We demonstrate the feasibil- 
ity and potential applications of this approach through 
3D finite difference time domain (FDTD) simulations, 
using the Lumerical 8 package. For the remainder of 
this paper, we assume that both the ring resonators and 
the access waveguides are made of 220nm high silicon 
(500nm width) encapsulated in Si02 and the two layers 
are separated by 275nm of silica. This layer separation 
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Figure 4. Simulated far-field intensity distribution of the 
beam emitted from a ring resonator, for (a) I = 1 and (b) 
1 = 2. The absence/existence of a doughnut pattern prove 
the vortex properties of these beam. However, the radial pat- 
terns show additional maxima due to angular diffraction for 
the finite sized source. 



was chosen as it provides good coupling to a resonance 
at 1542nm, our design wavelength. A top silica layer was 
included, to provide a comparable structure to those re- 
ported in Ref (27| . 

As a first step, we consider a simple case, where a single 
ring resonator couples to an access waveguide vertically, 
demonstrating that the vertical coupling does not affect 
the generation of OAM beams. This ring resonator has a 
radius of Rq = 3.9/im and the resonance at 1542nm has 
mode number m = 39. Two instances of this ring res- 
onator with gratings consisting of q\ = 38 and qi = 37 el- 
ements were simulated, giving an angular phase matching 
(and therefore OAM) of l\ = 1 and l<i = 2, respectively. 
Figure $4$ shows the simulated intensity patterns of the 
emitted beams in the far-field zone. The results con- 
firm that the emitted optical field for the case of £2 — 2 
possess a vortex structure, and in comparison to the hor- 
izontal coupling demonstrated in Ref [27]], the emitted 
optical field is not affected by the change to a vertically 
displaced access waveguide. For the case of l\ = 1 the 
emitted beam has a non-zero intensity at the centre, as 
discussed in the previous section. Apart from higher ra- 
dial mode excitations, the simulated far-field patterns are 
in excellent agreement with our theoretical model shown 
in Fig. ©. 

While a single ring resonator emits one of the complex 
polarization structures shown in Figs. @ and (U) this is 
not yet an OAM qudit state. In order to generate and 
manipulate OAM qudit states, we must address several 
concentric ring resonators at the same time. These ring 
resonators must be coupled to separate access waveguides 
and must have a resonance at the same wavelength, such 
that a superposition state can be formed. Figure ([5]) 
shows the resonance spectrum for three concentric ring 
resonators. The resonators are designed such that they 
all have a resonance at 1542 nm. When a single resonator 
is addressed, we observe a low crosstalk ( 10 dB) between 
the resonators. Due to the excessively large simulation 
time required for the three resonator system, most cal- 
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Figure 5. Spectrum of three concentric ring resonators: green 
dot-dashed, red dashed and blue solid lines correspond to in- 
ner, middle and outer resonators, respectively. All resonators 
have a resonance at 1542 nm, allowing for the formation of a 
superposition state. The resonator Q factors are between 700 
and 1000. 



dilations in this work are performed for a two resonator 
system. However, the results are general and can be ex- 
tended to systems with three or more ring resonators. 



APPLICATIONS 

To demonstrate the potential of the proposed system 
we will now discuss several potential and interesting ap- 
plications. 

(i) Quantum information: implementing a novel se- 
cure quantum key distribution (QKD) protocol in the 
OAM state space requires the availability of higher di- 
mensions. In this case different OAM bases of a sin- 
gle photon can be used as an individual alphabet to en- 
code the information [20l |. As shown previously, a ring 
resonators can be used to emit a variety of different 
OAM values along a common propagation axis, running 
through the ring's centre, providing access to the higher 
dimensions of the OAM space. When two or more con- 
centric resonators are simultaneously excited each will 
emit its own vector vortex beam, independent of the 
other resonators. Addressing these ring resonators co- 
herently leads to a coherent superposition of the emitted 
vector vortex beams, i.e. ~Y^t c i since the rings have 
matched resonances, assuming an identical input polar- 
ization. Therefore, without loss of generality, the polar- 
ization of the output beam can be factorized out. Let 
us now consider the simplest case, where the information 
is encoded onto two mutually exclusive bases of OAM. 
Such a bi-dimensional state is known as a qubit with 
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Figure 6. Simulated intensity pattern of a superposition state generated by two concentric ring resonators, (a) The beam is a 
superposition of I = +2 and £' — —2, and (b) I — +2 and £' = +1, respectively. In (b) the relative phase of the input signals 
for the two waveguides is varied, leading to a rotation of the superposition pattern. In (b), the intensity patterns, from left to 
right, correspond to relative phase of \ = 0°, % — 120°, and x — 240°, respectively. 



where (r\£) = e ltv is the OAM value, and 8 and x 
are the altitude and azimuth angles of the state in the 
Bloch sphere, which define the qubit's amplitude and a 
relative phase, respectively. Analogous to polarization 
states, the bi-dimensional OAM superposition state can 
be mapped over an OAM Poincare sphere, where the 
equator describes an equal superposition state having 
petal shapes with different relative phases [34| . 
To demonstrate that our system can be used to generate 
and control such states, we simulate two different 
superposition states. The first is a state consisting of 
equal magnitude yet opposite sign OAM values, more 
specifically £' = — £ = 2, while the second has different 
magnitude yet equal sign OAM values, £ = 2 and £' = 1. 
Figure [6]-(a) shows the first case, where two resonators 
are excited coherently and the gratings are designed 
such that the inner ring and outer ring emit a beam 
with £ = +2 and £' = —2, respectively. As can be seen 
from the far field pattern, the emitted beam forms a 
petal shape similar to a Hermite-Gauss mode of HGi.i 
(the emerging beam is a superposition of exp (2iip) 
and exp(— 2iip), i.e. cos 2^ with additional angular 
diffraction effects). In addition to forming superposition 
states, the relative phase of the two components of the 
superposition state, i.e. %, can be controlled. Figure 
[5]-(b) shows the simulated intensity pattern for the 
second case, where the inner and outer rings emit a 
beam with £ = +2, and £' = +1 respectively. For this 
beam the intensity pattern is known as a "C-beam". 
If a phase delay is introduced in the access waveguide 
of one resonator, then the relative phase between two 
OAM states in the superposition changes. Therefore, the 
C-beam intensity pattern rotates, as shown in Fig. |U-(b). 

(ii) Angular phased array: as shown in Fig. |6]-(b), 
the relative phase of the emitted beams can be adjusted 
through a control of the relative phase of the input 
waveguides. Interesting effects also occur if this idea 



is applied to concentric ring resonators emitting the 
same vector vortex beam, i.e. same wavelength, OAM 
value and polarization. Such rings form an angular 
array of coherent sources that are separated by a few 
wavelength, analogous to a phased antenna array. In 
this angular phased array, changes in the relative phase 
of the different emitters lead to an angular displacement 
of the emission pattern, demonstrated through a change 
in the relative intensity of the different radial orders. 

(iii) OAM sorter: quantum communication and QKD 
applications using OAM states require the ability to 
read out the encoded information and therefore an 
OAM sorter. In order to provide a scalable, i.e. mass 
producible system, an integrated, small scale sorter is 
required. All available sorters have a relatively large 
singularity size, on the order of a few micrometers, and 
require bulk optics components, which cannot be used 
during on chip integration (35| - |37| . However, due to 
the time (position) reversal invariance of the Maxwell's 
equations, a grating can act as both an emitter from and 
a coupler to a chip. Therefore, our resonator structures 
should not only be able to emit beams carrying OAM 
- they should also couple such beams into the access 
waveguide of the resonator. This effect was postulated 
for a single ring resonator in Ref [27j]. As shown in Fig. 
(|7|) , coupling does indeed occur if the rings are excited by 
a suitable vortex beam at vertical incidence, assuming 
that the optical axis is running through the center of the 
ring resonators. Furthermore, the coupling is selective, 
i.e. a ring resonator is only excited if the wavelength and 
the OAM value of the incident beam match the design 
of the ring resonator. As such, the OAM of an incoming 
beam can be measured by simply observing the output 
of ring resonators of known designs. If several concentric 
rings are used then a larger alphabet of values can be 
detected simultaneously and superposition states can 
be decomposed into their individual OAM components. 
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As a drawback, it is worth mentioning that there is a 
non-zero (— lOdB) crosstalk among different OAM values 
in the sorter configuration. A reduction of this crosstalk 
will require a design optimization, to be addressed in 
future work. 

CONCLUSIONS 



We present an approach to generate and manipulate 
the OAM of light in an on-chip setup, based on multiple 
concentric ring resonators. It has been shown that the 
emitted beam from the present configuration as well as 
other ring resonator based approaches form complex 
patterns of polarization, known as vector vortex beams. 
To our knowledge this setup is the first that allows the 
generation and dynamic control of OAM superposition 
states in an integrated fashion. Therefore, we believe 
that this design is a key component in the quest to 
encode and decode information using the OAM of light. 



Figure 7. Optical output power of the access waveguides for 
ring resonators that are set up for I — +1 (inner ring) and 
£ — +2 (outer ring) at 1542 nm. Red dashed and blue solid 
curves correspond to output power from the access waveguide 
of inner and outer ring resonators, respectively, (a) Shows the 
output beam power from access waveguides when illuminated 
with an equal superposition of I = +1 and £ = +2. The 
power unbalance is due to different out-coupling efficiencies 
for the outer and inner resonators. We can see a clear peak 
at the resonance of 1542nm. (b) Output power when illumi- 
nated with a superposition of £ — —2 and £ = +2 beams. 
The matched ring resonator (outer) couples at the correct 
resonance, while the other ring resonator (inner) couples at 
resonances for which the angular phase matching condition 
is met. Here we should note that the large bandwidth of 
the illumination is necessary for a shorter FDTD simulation 
time, however, experimentally the input could be chosen such 
that only the resonance of interest (here 1542nm) would be 
excited. 
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